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We report on the properties of a dipolar, antiferromagnetic chain in the framework of linear 
spin-wave theory. As a model we use an isotropic Heisenberg-Hamiltonian with antiferromagnetic 
nearest neighbor exchange and anisotropic dipole-dipole interaction. The phase diagram is calculated 
for an arbitrary ratio of dipolar interaction to exchange interaction and for fields both parallel and 
perpendicular to the chain direction. We can distinguish two regions, an antiferromagnetic exchange 
dominated region, in which the dipole energy acts as a planar anisotropy and a ferromagnetic dipole 
dominated region, in which the dipole energy acts as an easy-axis anisotropy. The anisotropy in 
the field dependent magnetization, which is caused by quantum fluctuations, is several orders of 
magnitude larger than the anisotropy predicted by the classical theory. We show that the dipole- 
dipole interaction does not lead to long-range order in the pure one-dimensional case, however for 
quasi one-dimensional systems spin- waves are important. The calculated magnetization is compared 
with experimental measurements for RbMnBr3, a quasi one-dimensional antiferromagnet. 
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I. INTRODUCTION 

Low-dimensional magnetic systems have long attracted 
the attention of experimental and theoretical physics. In 
systems with reduced dimensionality the fluctuations are 
enhanced and can prevent a system from long-range or- 
dering. For one- and two-dimensional isotropic Heisen- 
berg ferro- and antiferromagnets,with short-range inter- 
action it has been proven exactlytra that long-range order 
is absent for finite temperatures. The existence of long- 
range order in low-dimensional systems depends crucially 
on anisotropies and long-range interactions, as there are 
single-ion anisotropy or the dipole-dipole interaction. 

In real systems the dipole-dipole interaction is always 
present in addition to the short-range exchange interac- 
tion. Although the dipole-dipole energy is weak com- 
pared to the exchange energy it plays an important role 
in low-dimensional systems due to its anisotropic and 
long-range character. MaleevB showed that in a two- 
dimensional ferromagnetic system the additional dipolar 
interaction leads to a modification of the spin-wave fre- 
quency in the long wavelength limit from a quadratic to 
a square root wave vector dependence, for which long- 
range order is established. For the antiferromagnetic 
counterpart on a square and a honeycomb lattice it has 
also been shownQ that dipolar interaction stabilizes long- 
range magnetic order at finite temperatures. In this 
cases, however, the dipolar anisotropy leads to a discrete 
ground state with spins oriented perpendicular to the 
plane, producing an energy gap in the dispersion rela- 
tion. The finite Neel temperatures in various quasi two- 
dimensional systems (e.g. .K2M11F4) have been calculated 
within the dipolar modclH, in good agreement with the 
experimental values. 

In recent papersB we determined the magnetic phase 
diagram of the two-dimensional, dipolar Heisenberg an- 
tiferromagnet (square and honeycomb lattice) for fields 
along the easy-axis, i.e. perpendicular to the plane. 
In contrast to models with singlc-ion anisotropy or an 
anisotropic exchange interaction, we found that the 
dipole-dipole interaction establishes a new phase, an in- 
termediate phase. With increasing field strength the sys- 
tem changes by second-order transitions from the Neel 
phase to an intermediate phase, a spin-flop phase and 
finally to a paramagnetic phase at high magnetic fields. 

From the results in two dimensions we anticipate that 
also in quasi one-dimensional Heisenberg antiferromag- 
nets, the dipolar interaction will reduce thermal and 
quantum fluctuations and thus increase the tendency to 
order. 

In this paper we investigate the influence of the dipolar 
interaction on the classical ground states and on quan- 
tum fluctuations. It is shown that due to the competition 
of the exchange and dipolar interaction we can distin- 
guish two regions, namely an exchange dominated and 
a dipolar dominated region. In the former the dipolar 
energy leads to a planar anisotropy, whereas in the lat- 



ter it leads to an easy-axis anisotropy. We study the 
dispersion relations for various relative strengths of the 
dipolar interaction. In the exchange dominated region, 
the dipolar interaction breaks the continuous symmetry 
of the ground state partially and favors states with spins 
perpendicular to the chain direction. Rotation symme- 
try around the chain direction still exists. The dispersion 
relation for this phase is modified by the dipolar inter- 
action, but as in the isotropic counterpart no long-range 
order can exist even for T = 0. 

The most famous quasi one-dimensional systems are 
the ternary compounds ABX3 (A alkaline, B transition 
metal and X halogen), which have been studied inten-, 
sively experimentally in ±he context of Haldane's phasecl 
and solitonic excitationsQ. Of particular interest are the 
Mn compounds. Because the angular momentum L is 
zero, no crystal field splitting occurs in these systems 
and the dipole-dipole interaction-should be the most im- 
portant anisotropy. Experimentaj show an anisotropy in 
the magnetization measurements for fields applied paral- 
lel and perpendicular to the spin-chain axis. We will show 
that this anisotropy can be understood qualitatively by 
our simple one-dimensional dipolar model. Up to now, 
only partial asp«jte of the dipole-dipole interaction have 
been considerecot2l. In Rcf[9] the classical ground state 
has been determined, but no excitations have been inves- 
tigated, whereas in Ref[10] the dipolar interaction is only 
considered for nearest neighbors. In this paper the effects 
of the additional dipolar interaction will be studied in a 
systematic way. 

In this spin wave theory we disregard kink-like excita- 
tions, which even for the r -3 decay of the dipolar inter- 
action destroy one-dimensional long-range order. How- 
ever, for wavelengths shorter than the correlation length 
magnon excitations become important. Tpbis has been 
shown experimentally by Steiner for CsNiFatU. The short 
range antiferromagnetic order supports spin waves with 
wave vectors not too close to the center of the Brillouin 
zone. Especially for quasi one-dimensional systems the 
weak interchain interaction leads to a stabilization of the 
magnetic order and spin waves become the relevant ex- 
citations. 

The outline of this paper is as follows: in section II 
we introduce the dipolar model; in III we investigate the 
ground state, the spin-wave spectrum and the order pa- 
rameter for vanishing fields; in IV we study the magnetic 
phase diagram for the whole parameter region. In the last 
section we apply our model to a quasi one-dimensional 
system, namely RbMnBr3. 

II. MODEL 

The Hamiltonian of a dipolar antiferromagnet reads 
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with spins S; at lattice sites x;. The first term in brackets 
is the isotropic exchange interaction Ju> and the second 
term is the classical dipole-dipole-interaction 
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The third term in Eq. (2.1) describes the effect of an 
external magnetic field Ho, where g is the Lande- factor 
and [Lb is the Bohr magneton. From this Hamiltonian 
we calculate the magnon spectoim in the framework 
of the Holstein-Primakoff (HP)E3 transformation. The 
Holstein-Primakoff transformation is an exact transfor- 



mation which expresses the Hamiltonian of Eq. (2.1) in 
terms of Bose operators a, , a; . For low temperatures an 
expansion up to bilinear terms can be used 



a ai 



(2.3) 



The tilde indicates that we have to use these inner spin 
coordinates in the rotated coordinate frame where the z- 
direction lies in the direction of the classical spins. For 
the following discussion we define a dimcnsionless param- 
eter k, which gives the ratio between the dipolar energy 
and the exchange energy: 
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Here J denotes the strength of the exchange energy (the 
negative sign of the antiferromagnetic exchange interac- 
tion Jui is taken into account separately) and c denotes 
the distance between magnetic ions. 
In the following we consider a spin chain which is directed 
along the x-axis (see Fig. |l^,b) . 
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where J q is the Fourier transform of the exchange inter- 
action Jui, which equals J q = — 2 J cos(q x c) for nearest- 
neighbor interaction (in the following we set q x = q), 
and is the Fourier transform of the dipolar tensor 



Eq. Q2.2D 
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which can be evaluated by Ewald summationliij (see ap- 
pendix A). The wave vector q = -(1, 0, 0) describes the 
antiferromagnetic order. Eq. leads to a critical 

value of k 
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with C(^) an d Tj{n) defined in appendix A. 

In the following we investigate the classical ground 
states and excitation spectra of the one-dimensional 
chain as a function of k. 



A. Exchange dominated region (k < k c ) 

For pure exchange energy J, the classical ground state 
is the configuration where the spins are antiparallel, as 
described by the antiferromagnetic wave vector qo . This 
state is continuously degenerate due to the rotational 
symmetry. The dipolar interaction breaks the rotational 
symmetry and favors the antiferromagnetic ground state, 
in which the spins lie in a plane perpendicular to the 
chain axis, as can be seen by the largest value of the 
dipole-dipole tensor (A^ o = A^ Q > A^ ). However, rota- 
tion symmetry around the x-axis still exists (s.Fig. 0a). 



(a) 



III. DIPOLAR ANTIFERROMAGNETIC CHAIN 
WITHOUT FIELD 

The long-range dipole-dipole interaction and the an- 
tiferromagnetic exchange interaction compete with each 
other: while the latter favors an antiferromagnetic spin 
configuration, the former favors a ferromagnetic orienta- 
tion along the chain axis. Thus, we can distinguish an 
exchange- dominated region and a dipole-dominated re- 
gion, w hich a re separated by the condition (to be shown 
in Eq. (Qgh) 



(b) 



FIG. 1. a) Classical ground state of the dipolar, antiferro- 
magnetic chain in the exchange dominated region (k < k c ). 
b) Ground state for the dipolar dominated region (k > n c ). 

In the following calculations the primitive cell is the 
chemical cell which is half the magnetic unit cell. The 
antiferromagnetic modulation is taken into account via 
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where £1(^2) denotes the first (second) sublatticeQ. If 
we start from the classical ground state with the spins 
oriented along the ±z-axis and perform the Holstein- 
Primakoff transformation, we get the following Fourier 
transformed Hamiltonian 

H = E c i + ^4 q a q a q + i_B q (a q a_ q + a q aL q ) (3.5) 



with the classical ground state energy 
S c/ = -iV5 2 (J qo +4 3 o ) 
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and the coefficients 
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The Hamiltoniaa (3__5) is diagonalized by a Bogoliubov 
transformation!!^ and reads in terms of the new Bose cre- 
ation and annihilation operators & q and 6 q 

H = E « + ^E^-^ + E^ 6 ^- ( 3 - 8 ) 



The second term in Eq. ( 3.8 ) is a correction to the ground 
state energy due to quantum fluctuations, while the third 
term is the harmonic magnon Hamiltonian. The ground 
state energy 

£(T = 0)=£ c , + ^(£: q -Aq) (3.9) 
q 

is shown in Fig. ^| as a function of k. The quantum fluc- 
tuations are strongest for pure exchange interactionll3. 
The deviation of the ground state energy from the clas- 
sical value 
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is most significant in this case and decreases with increas- 
ing dipolar interaction, i.e. the dipolar interaction tends 
to stabilize the antiferromagnetic ground state. 
The magnon dispersion relation reads 
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which is shown in Fig. || for four values of k. Because of 
the continuous degeneracy of the ground state (rotation 
around the x-axis) the Goldstone-mode at q = sur- 
vives in the presence of the dipolar interaction, but the 
behavior for small wave vectors is modified. If we expand 
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FIG. 2. Deviation of the ground state energy from the clas- 
sical value due to quantum fluctuations for k < k c - For k > ti c 
there is no difference between the two energies, because the 
Hamilto nian is already diagonal after the HP-transformation 
(s. Eq. (E20k). 
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FIG. 3. Dispersion relation for different values of k in the 
exchange dominated region (k < k c ). 



the magnon spectrum around the center of the Brillouin 
zone, we get logarithmic corrections 

Ej(iJS) = U(qc) 2 - + {qcf log (qc) 



+b 2 (qc) i + 0((qc) 



1/2 



(3.12) 



with g-independent constants C5, b%, 62 (see Eq. A7, A14, 
A15 in Appendix A). 



The order parameter, the staggered magnetization, is 
given by 
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with 07 from Eq. (3.4). Inserting the Holstein-Primakoff 
transformation (2.3), we obtain 
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with the zero-point contribution 
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and 
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from thermal excitations of spin waves. n q = Eq /i T — - 
denotes the Bose-occupation number. Both quantum and 
thermal fluctuations reduce the order para meter . By use 
of the small wave vector expansion, Eq. ( 3.12| ), we see 
that the sum No diverges for q — » 0. Even though the 
dipolar interaction has the tendency to stabilize the anti- 
ferromagnetic order, it is nevertheless destroyed by quan- 
tum fluctuations even at T = 0, as found for the isotropic 
counterpart. There is yet another stability condition re- 
sulting fro m th e requirement that the spin wave frequen- 
cies (Eq. ( |3 .11 )) be real which leads to 
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for all wave vectors We get a lower bound for the 
stren gth of the exchange energy from the first condition 
(Eq. at q = 
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For lower exchange energies the system changes to the 
dipolar dominated region. 



B. Dipolar dominated region (k > k c ) 

For pure dipolar interaction the classical ground state 
consists of two states, in which the spins are oriented 
ferromagnetically along the ±x-axis (Fig. |l|b), because 
the largest eigenvalue of the dipolar tensor is Aq 1 (s. 
Fig. i.e. the dipolar interaction acts as an easy- 

axis anisotropy. If we start with this ground state and 
perform the HP-transformation we get 



with the classical ground state energy 
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Due to the broken rotation symmetry the spectrum has 
an energy gap at the zone center 
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FIG. 4. Spin-wave dispersion relation of the dipolar, 
tiferromagnetic chain for k — k c = 2.66 (solid) and K 
(dashed). 
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The minimum of the dispersion relation is found at the 
zone edge (s. Fig. 4), which vanishes for n — n c and 
increases with increasing dipole strength. 

In linear spin-wave theory the classical ground state 
and the quant um m echanical ground state coincide. The 
Hamiltonian ( 3.20|) is diagonalized by the HP-trans- 
formation, .and there are no corrections for quantum 
fluctuationsEEl at T = 0, i.e. the ground state has long- 
range order. The order parameter is the magnetization 
in x-direction: 

M x = g^B ^{Sf) = g^sNS - M th (T) (3.23) 



with 
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For k > k c the spin- wave spectrum has gaps at q = and 
q = q (Eq. Q3.22fl). Thus, the sum in Eq. ( p\2l| ) is finite 



for finite temperatures and from point of view of the lin- 
ear spin- wave theory long-range oider should be possible. 
However, as was shown by DysonLj, kink-like thermal ex- 
citations destroy the long-range order in ferromagnetic 
Ising systems with interactions decreasing like J(r) oc ^ 
for a > 2. Thus, we expect that for a Heisenberg system 
with dipolc-dipole interaction these nonlinear excitations 
will ultimately destroy the long-range order as well. 

Stability of the ground state requires that the spectrum 
is positive, Eq > 0. For decreasing dipole energy the 



instability of the spectrum sets in at an, 
recover the stability condition Eq. ( |3.l[ ) 
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The soft mode at qo for k = n c indicates the phase transi- 
tion to the antiferromagnetic phase (exchange-dominated 
region) . 
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Consequently at T = there is long-range order (LRO) 
in the dipolar dominated region, in contrast to the ex- 
change dominated region. However, for finite tempera- 
tures no LRO exists at all due to nonlinear excitations, 
namely domain walls. 

IV. MAGNETIC PHASE DIAGRAM 

In this section we determine the magnetic phase dia- 
gram of the dipolar antiferromagnetic chain for both the 
dipolar dominated and the exchange dominated region. 
We study homogeneous magnetic fields along and per- 
pendicular to the chain axis (^-direction). The r esul ts 
are presented in terms of the parameter k (Eq. (EM)) 
and the dimensionless field h — gp f s " . The spin ori- 
entations are determined from the classical ground state 
energy within a two-sublattice model. 




FIG. 5. Magnetic phase diagram of the dipolar, antiferro- 
magnetic chain with a field applied in chain direction. The 
spin arrangements below the a;-axis indicate the phases for 
vanishing field. In the insert the configuration of the spins in 
a two-sublattice model is shown. 



A. Field in chain direction 

First we study the case of a magnetic field along the 
chain axis. In the dipolar region (k > k c ) the field is par- 
allel to the spins, and thus no reorientation takes place. 

In the exchange dominated region (k < k c ) the spins 
are oriented perpendicular to the x-axis and any infinites- 
imal field along the axis will lead to a reorientation of the 
s pin s in order to gain energy from the Zeeman term (Eq. 
( [2.l|) ). In a two-sublattice model (the spins of the sublat- 
tice £i(£2_) enclose an angle <^i(<y9 2 ) with the y-z-plane 
(see Fig. [5])) we can calculate the classical ground state 
energy. For spins lying in the sz-plane we get 



The complete phase diagram is shown in Fig. ||. The 
HP-transformed Ham iltonian in the spin-flop phase has 
the form of Eq. fl3.5|) with coefficients 
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where tp is related to the field by Eq. (4.2). The dis- 
persion relation is shown in Fig. |6| for two values of the 
field. 
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Minimization with respect to ipi and <p 2 yields tp\ = <p2 = 
ip, i.e. a usual spin-flop phase. From the minimum one 
determines the relation between the angle ip and the mag- 
netic field: 
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Thus the spins become progressively aligned along the 
field as the field increases. There is a continuous transi- 
tion from this spin-flop phase to the paramagnetic phase 
for p = 90°. The transition occurs at 
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FIG. 6. Dispersion relation for the linear chain in the 
spin-flop phase for two values of k; solid lines for vanishing 
field and dashed lines for simp = 0.5 (Eq. (4.2)). 



The Goldstone-mode at q = remains for finite dipo- 
lar energy because the system is still invariant for rota- 
tions around the x-axis. 

The magnetization in the spin-flop phase can be cal- 
culated from the free energy F 
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The free energy F is related to the partition function Z 
by 



F=-kT log Z with Z = Spe- /3H and (3 
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i.e. for the Hamiltonian ( |3.8] ) 
F = ^ + | E^i - ^ + ^T^log(l - e-^) 
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Inserting Eq. (h6) into (4J3) yields for the magnetization 



M = M c; -M -M th (T) , 
which consists of a classical part 

dE cl 



M r 



9H 



a part arising from quantum fluctuations (T = 0) 
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and a part arising from thermal fluctuations 
1 1 



M th (T) = J2- 



■P E « - 1 E a 



(4.10) 



The classical magnetization is reduced by quantum and 
thermal fluctuations. For the spin-flop phase the classical 
magnetization can be obtained by inserting the expres- 
sion for the classical ground state energy (Eq. (|4.l|)) in 
Eq. (fy|). We get 



M% = g^FsNS sin ip , 



(4.11) 

where siny is given by Eq. ( |4.2| ). The contributions 
of quantum and thermal fluctuations will be taken into 
account numerically when we compare our results with 
experiments (s. section IV). 

Now we turn to the staggered magnetization in the 
z-direction 



N 2 

which leads to 
N z = g^B cos ip 
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The second and third term in (4.13) are divergent due to 
the Goldstone-mode. Thus, no antiferromagnetic long- 
range order is possible in this case, only a finite magne- 
tization. 



B. Field perpendicular to the chain direction 

1. Exchange- dominated region (k < k c ) 

A field perpendicular to the chain leads to a spin-flop- 
transition; for infinitesimal fields the spins lie in the yz- 
plane and perpendicular to the applied field. With in- 
creasing field the spins become progressively aligned in 
field direction. In the following calculations the field is 
applied in the y-direction so that for infinitesimal field 
the spins are oriented in the indirection. Again, we 
study the classical ground state energy of a general two- 
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FIG. 7. Configuration of the spins for a field applied per- 
pendicular to the chain axis; a) for the exchange-dominated 
region, k < k c , b) for k > n c . 

sublattice system, for which the usual spin-flop phase 
turns out to be stable. The equilibrium condition leads 
to 



sin ip 



2S{Jq -J Q +A™-Af) 



(4.14) 



in analogy to Eq. (4.2) for fields along the chain axis, 
which differs only in the dependence on the dipole inter- 
action. The transition between the spin-flop phase and 
the paramagnetic phase occurs at (see Fig. [l(]) 



h v c {K) =8 + k(2?7(3) + 2C(3)) for 



(4.15) 



The coefficients of the Hamiltonian, Eq. (3.5), in the 
spin-flop phase are in linear spin-wave theory 
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For finite k rotational symmetry is lost for the spin config- 
uration of Fig. ^a and thus no Goldstone-mode is present 
in the dispersion relation of Fig. ^ showing the spectrum 
for two different values of k and H^. The magnetiza- 
tion in the y-direction is obtained from the free energy 
(see Eq. (4.5)) and is of the same form as Eqn. fl4.7| )- 
( 4.11 ) with different canting angle sin t p (Eq . (| 4. 14| ) ) and 



different coefficients Aq and Bq (Eq. (4.16)). 
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h v c {n) = 6£(3)k for k > k c 



(4.20) 
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The coefficients in the Fourier transformed Hamiltonian 
for finite canting angle d are given by 



A q = 2S(J - J q ) + 2S(A 1 1 - A 22 ) 
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The spectrum is shown for three values of the field in Fig. 

For decreasing field the energy gap Eq decreases until 
it vanishes at the transition point. Now we investigate 
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FIG. 8. Dispersion relation in the spin-flop phase for differ- 
ent magnetic fields applied perpendicular to the chain direc- 
tion and different strengths of the dipolar interaction. The 
solid lines are w itho ut field and the dashed lines are for 
simp = 0.5 (Eq. (O). 



Now the staggered magnetization in z-direction, which 
in analogy to Eqs. ( |4. 1 2| )- ( p. 13| ) is given by 
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is finite, because of the energy gap in Eq (s. Fig. g). 
However, we believe that due to nonlinear kink-like ex- 
citations again no long-range order is possible for finite 
temperatures. 

If the field is applied in chain direction a Goldstone- 
mode remains (s. Fig. [|) and hence the expression for 
the staggered magnetization N z diverges (Eq. fl4.13| )). 



2. Dipole dominated region (k > n c ) 

In this region the spins keep the ferromagnetic orien- 
tation but with a canting angle $ (intermediate phase) 
with respect to the magnetic field (s. Fig. 7b). The 
classical ground state energy is given by 
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-g(i B H%NS sintf , 



i? + ( Jo + Al 1 ) 



cos 



(4.18) 



the minimization of which gives the following relation 
between the canting angle -d and the applied field Hq : 
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Note that this angle depends only on the dipole strength 
but not on the exchange energy due to the isotropy of 
J q . The transition to the paramagnetic phase, in which 
all spins orient along the magnetic field, occurs at 




FIG. 9. Dispersion relation for the dipole dominated re- 
gion with a field perpendicular to the chain direction and 
k = 4; dashed: sini? = 0, dot-dashed: sun? = 0.5 and solid: 
sin i? = 1. 

the transition from the spin-flop phase to the intermedi- 
ate phase as a function of k. Comparison of the classical 
ground state energy of the two phases shows that the 
transition line is independent of the applied field (if 
Hq < H%), which explains the vertical line at k = k c 
between the spin-flop phase and the intermediate phase 
in the phase diagram of Fig. 



10. For k < k c the spin- 



flop phase has lower energy and is thus the stable phase, 
for k > k c the intermediate phase has lower energy. At 
k = k c and h = h c the energies of the spin-flop phase, 
the intermediate phase and the paramagnetic phase are 
the same. 



V. COMPARISON WITH EXPERIMENT 

In this section we compare our theoretical results with 
measurements on quasi one-dimensional antiferromag- 
nets in the hexagonal ABX3 structure. Due to different 
exchange mechanisms the intrachain interaction is about 
three orders of magnitude larger than the interchain in- 
teraction, thus that these materials behave in a quasi one- 
dimensional manner. The most promising substances 
are the manganese compounds, because Mn 2+ ions have 
spin S = 5/2 and vanishing orbital momentum. Thus, 
the main contribution to the anisotropy should be the 
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h^(K c ) 



h c y (0) 



FIG. 10. Magnetic phase diagram of the dipolar, antiferro- 
magnetic chain with a field perpendicular to the chain axis. 
The configurations below the z-axis correspond to vanishing 
field. 



dipole-dipole interaction. The field-dependent magneti- 
zation, has been measured for RbMnBr3 (Fig. |ll| dashed 
lines)Q. In this material the spins are forced to orient 
within the hexagonal plane, perpendicular to the chain- 
axis. This planar anisotropy can be fully explained by the 
dipole-dipole interaction, which leads in the exchange- 
dominated region to just such an orientation (see section 
II). Magnetization experiments show an anisotropy: the 
magnetization for magnetic fields along the spin-chain 
axis is about 5 to 10% larger than for fields perpendicular 
to it. Part of this anisotropy can be understood qualita- 
tively even for classical spins by the dipolar anisotropy, 
which leads to different dependenc es of the a ngle ip of the 
spins on the magnetic field (Eqs. (4.2) and (4.14)). The 
classical contribution to the magnetization is larger for 
parallel than perpendicular field, since the dipole interac- 
tion favors ferromagnetic orientation in chain direction. 
For a quantitative comparison we need the values of the 
exchange energy J and-the lattice constant c which are 
given for RbMnBr3 byE] 



J = 9.56^, 



3.26A . 



Here the exchange energy was determined by a fit of the 
spin- wave curves with the spin-wave spectrum calculated 
from a Hamiltonian in which the anisotropy is described 
by a single-ion term. For this Hamiltonian there exists a 
renormalization calculationEJ, which shows, that the clas- 
sical value of the exchange energy has to be renormalized 
due to magnon-magnon interactions. The resulting value 
of J is 

J = 8.93^ . 

This corresponds to k = 8.07 • 10~ 3 , which is deep in 
the exchange dominated region of our theory. With this 
constant we calculate the class ical contributions to the 
magnetization (see Eq. (4.11)) in the spin- flop phase 
shown in Fig. [ll] as the upper and lower thin solid 



lines for fields along and perpendicular to the chain- 
axis. The dashed lines represent the experimental values. 
The transition to the paramagnetic phase takes place at 
H c k 130T and cannot be seen here. For perpendic- 
ular fields lower than H « 3.9T the three-dimensional 
ordering becomes important: the system changes to a 
complicated six-sublattice structured for which our one- 
dimensional model is inappropriate. According to Fig. 
[Til , the classical magnetization is larger than the experi- 
mental and the splitting for the two directions is almost 
invisible. However, in one dimension the quantum fluc- 
tuations are strong and must be taken into account. Nu- 
merical evaluation of the magnetization for longitudinal 
and perpendicular field direction (Eq. (|4.7|)), including 
the zero-point contributions (T = 0), yields the thick 
solid lines in Fig. [H]. In fact, the quantum fluctuations 
lead to a pronounced reduction of the magnetization and 
increase the splitting, but as to be expected, are some- 
what overestimated in the pure one-dimensional model. 
The contribu tion of thermal fluctuations to the magneti- 
zation (Eq. ( [4.10 )) is negligible for temperatures as low 
as in the experiment (T p,1.7A"). This is in contrast to 
the work of Santini et al.ll3, who claim that the classical 
thermal fluctuations are the relevant contribution and 
totally nfijglect the quantum fluctuations. Abanov and 
PetrenkoQ also studied the one-dimensional spin-chain, 
but with a single-ion anisotropy favoring orientation per- 
pendicular to the chain instead of the dipolar interaction. 
Inserting the fitted value for this anisotropy they get a 
larger fluctuation contribution than in our theory. Fur- 
thermore, the classical contribution shows the opposite 
anisotropy, i.e. the magnetization for a field along the 
chain-axis is lower than for a perpendicular field. 

In summary, we can explain the anisotropy measured 
in RbMnBr3 qualitatively by a one-dimensional, dipolar 
Heisenberg model. We expect that by taking into account 
the full three-dimensional structure with interactions be- 
tween different chains the quantitative agreement could 
be improved because the contribution of quantum fluc- 
tuations is reduced in higher dimensions. Our prediction 
for the spin wave dispersion relation of RbMnBr3 (Eq. 

P) is shown in the magnetic Brillouin zone in Fig. 
[ere we set the exchange interaction J to 9.56 K. 
The resulting gap of the upper dispersion branch at q = 
equals 14.16 K. 



VI. SUMMARY 

We have considered the one-dimensional, isotropic 
Heisenberg antiferromagnet with dipole-dipole interac- 
tion in the complete parameter region. We showed that 
there exists an exchange dominated phase for k < k c , 
where the spins orient antiferromagnetically, perpendic- 
ular to the chain axis, i.e. the dipole-dipole interaction 
gives rise to a planar anisotropy. For k > k c , the dipole- 
dominated phase, the ground state is ferromagnetic along 



9 



experiment 

classical theory 

theory with quantum fluctuations 




2.0 3.0 

field (tesla) 



7I/2C 



FIG. 12. The two branches of the dispersion relation for 
RbMnBr3 in the magnetic Brillouin zone with dipolar inter- 
action (k = 8.072 ■ 1(T 3 , J = 9.56 K). 



FIG. 11. Magnetization curves for RbMnBr30: The dashed 
lines show the experimental measurements at T = 1.7K. The 
thin lines show the classical results of our theory, i.e. without 
quantum fluctuations. The thick lines show the results taking 
into consideration quantum fluctuations. For each pair of 
curves the upper one refers to field parallel and the lower one 
to field perpendicular to the chain. 



the chain axis, i.e. the dipole-dipole interaction acts like 
an easy-axis anisotropy. Because of the broken rotation 
symmetry the dipole-dominated region shows long-range 
order for T = 0. In the exchange-dominated region there 
is still a Goldstone-mode due to the rotation symmetry 
around the chain axis. Although the linear small q be- 
havior of the acoust ic ma gnon mode is modified by the 
dipole energy (Eq. ( 3.12) )), there is no long-range order 
neither for finite nor for zero temperature. 

The magnetic phase diagram for fields along and per- 
pendicular to the chain in both the dipole and the ex- 
change dominated regions has been determined and the 
predicted magnetization has been compared with exper- 
iments on RbMnBr3, a quasi one-dimensional, planar 
antifcrromagnet. The experimental anisotropy can be 
explained qualitatively by the one-dimensional, dipolar 
model, the quantitative agreement could be improved by 
a more elaborate, three-dimensional calculation. 
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APPENDIX A: DIPOLAR SUMS 

The dipolar surns.can be evaluated in Fourier space by 
the Ewald methodEl. The sum in Fourier space is divided 
into a sum over the direct lattice and a sum over the 
indirect lattice, so that both sums converge quickly. The 
dipolar tensor in one dimension finally reads (all values 

of the dipolar tensor are measured in units of 



) 



A<*P 



|7r-27rE ; ^ x Vi(^|x ; | 2 ) 
-27rE r ^- 2 (^|q+G r | 2 ) a = /3 # x, 

|tt5 q/3 - 2nS aP j:' l e'^Vt (£|x,| 2 ) 
+4^Eie 4qX! ^(V i (£|x ; | 2 ) 
-c 2 E r (q+G r ) a (q + G r )^ 

a, (3 £ x/ , 
(Al) 



where x, = Zc(l,0,0) and G r = ^r(l,0,0). The prime 
at the /-summations denotes I ^ 0. We have to dis- 
tinguish between components with spins lying in the 
chain direction (a, (3 6 Xj) and those perpendicular to 
it (a, (3 £ Xf). The off-diagonal components are iden- 
tical zero due to the inversion symmetry. The dipolar 
sums are shown in Fig. (|ll|) as a function of q. In Eq. 
(Al) we have introduced the Misra functions 



dtt n e~ xt 



(A2) 



The asymptotic behavior of the Misra functions for small 
arguments x is 



<p-i(x) w (-7 -logs) +. 



1 



0{x 3 



(A3) 



and 
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ip- 2 {x) as 1 + (-1 + 7 + loga;)x - -x 2 + 0{x 3 ) , (A4) 

where 7 « 0.5772 is Euler's const ant. The expansion of 
the dipolar tensor needed in Eq. ( 3.12 ) reads 



.4, 



q - Ci - c 2 {qc) + (qc) log 
+c 3 (qcf + O [(qc) 6 ) 



47T 



(A5) 



4+qo = = c 4 - c 5 ( 9 c) 2 + c 6 (qc) 4 + O [(qcf) , 

(A6) 

where the cons tants c\ , C4 are given analytically in Eq. 
(AS) and (A10), while the others are computed numeri- 
cally as 

c 2 w 0.666 , c 3 « 0.066 , c 5 w 0.693 and c 6 « 0.021 . 

(A7) 



The dipolar tensor can be evaluated exactly for some 
configurations. For the ferromagnetic order of the spins 
along the x-axis one gets 



<= 4 E^=4C(3) 



d as A, 



(A8) 



where Q(n) is the Riemann zeta function0 

00 

C(n)=Er n . (A9) 



The antiferromagnetic order of the spins along the 2-axis 
yields 



4 3 = 2E 



(-1) 



1-1 



= 277(3) = c 4 w 1.803 (A10) 



where the eta function r)(n) is defined by0 

00 

^^(-r/- 1 /-. (All) 



1=1 



Other exact values are 

4f = -2C(3) and ^ = -4f,(3) . (A12) 



These values corroborate the sums in Eq. (Al). One can 
show that for the dipolar tensor components of the linear 
chain the relation 



is valid. 

The constants in Eq. ( |3.12 ) are given by 



(A13) 




FIG. 13. Dipole sums A^f in Fourier space given in units of 

(9 H B f . The solid line denotes Al 1 and the dashed A 2 a 2 
All non-diagonal components vanish. 



l(8-«(4C(3)-2t ? (3)))(2 + «c B ) 



+ K 



(A14) 



^(4C(3)-2r ? (3)) + i( C2 ~ C5 )-i C6 



1 

16 
1 



KC2C5 



1 

16 



Kc 6 (4C(3)-2r ? (3)) 



log(47r) 



1 

16 



kc 5 log(47r) 



(A15) 
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